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Abstract 

The subject of moving curves (and surfaces) in three dimensional 
space (3-D) is a fascinating topic not only because it represents typi- 
cal nonlinear dynamical systems in classical mechanics, but also finds 
important applications in a variety of physical problems in different 
disciplines. Making use of the underlying geometry, one can very 
often relate the associated evolution equations to many interesting 
nonlinear evolution equations, including soliton possessing nonlinear 
dynamical systems. Typical examples include dynamics of filament 
vortices in ordinary and superfluids, spin systems, phases in classi- 
cal optics, various systems encountered in physics of soft matter, etc. 
Such interrelations between geometric evolution and physical systems 
have yielded considerable insight into the underlying dynamics. We 
present a succinct tutorial analysis of these developments in this arti- 
cle, and indicate further directions. We also point out how evolution 
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equations for moving surfaces are often intimately related to soliton 
equations in higher dimensions. 

Running Title: Nonlinear Dynamics of moving curves and surfaces 
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1 Introduction 



Curves in 3-D correspond to one of the simplest set of objects imaginable in 
geometry. Commonly encountered physical systems which can be modeled 
on curves vary over a spectrum of subjects ranging from physical, chemical 
to life sciences. Examples include vortex filaments in ordinary and 
superfluids [Hasimoto, 1972; Schwarz, 1998], scroll waves over fluid surfaces 
[Cross and Hohenberg, 1993], polymers [Bird et al., 1977] and biological 
macromolecules such as proteins and DNA. A not too obvious example is 
the classical version of Heisenberg ferromagnetic spin lattice in its 
continuum limit [Lakshmanan et al., 1976], which will be discussed in detail 
later. All these above stated examples pose challenging problems in physics 
with regard to their dynamics due to their highly nonlinear character. 
Besides these examples, the curve as an abstract geometrical object poses 
several interesting questions associated with its geometry for the 
mathematically inclined. The subject has held the attention of physicists, 
applied mathematicians and the nonlinear dynamics community in 
particular, for, as we shall see, the equations governing their dynamics are 
inherently nonlinear. In certain physical systems the models closely 
resemble soliton possessing integrable systems. These models prove to be a 
useful starting point in unraveling their apparently complex behavior. 
Further, the ideas developed for curves provide a good lead to their 
generalization and extension to higher dimensional geometric objects such 
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as surfaces and hyper-surfaces. 

In this article, after a brief description of the geometry of moving curves, 
we review some of the interesting intrinsic quantities associated with them 
that are frequently encountered, followed by a pedagogical analysis of 
several applications to physical systems wherein substantial results have 
been obtained. We also indicate possible generalizations to surfaces. It is 
the aim of this article to point out that it is possible to study the dynamics 
of several physical systems with inherent nonlinearity at a fundamental 
level by identifying their underlying geometry. 
The paper is organized as follows. In section 2 we discuss the basic 
geometry of a curve in 3-D and introduce the Frenet-Serret equations which 
describe it. In section 3 we further discuss a few (global) geometrical 
quantities that are frequently encountered. Connection to moving curves in 
fluid mechanics has yielded significant results in the past. This will be 
discussed in section 4. Section 5 discusses the dynamics of the 1-D classical 
Heisenberg ferromagnetic spin chain in its continuum limit which has been 
shown to be integrable through its connection with curve motion. Polarized 
beams propagating along a curved path acquire a geometric phase 
depending on the geometry of the path taken. This phase is discussed in 
section 6. Recently there has been several results in the physics of soft 
matter that can be modeled on curve motion, which will be discussed in 
section 7. Section 8 discusses the extension of these ideas to the geometry 
of surfaces. We conclude in section 9 with a summary and discussion. 
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2 Geometry of Moving Curves and Soliton 
Equations 

Consider a curve in 3-D represented by the position vector R(#), where 9 is 
some parameter along the length of the curve. It is possible in principle to 
switch over to the arc length parameter s, such that the tangent vector 
ei(s) defined as (see Fig. 1) 

e,( S ),f (1) 

is of unit length all along the curve. With this definition of the unit tangent 
vector one can canonically define a unit normal vector e 2 (s), orthogonal to 
e 1; such that 

ei s = ree 2 . (2) 

Here (and henceforth), the subscript denotes derivative with respect to s. 
Notice that the derivative of e x does not have a component along ei itself 
owing to its constant modulus. This follows for the unit normal e 2 too, and 
also for the unit binormal e 3 to be defined next. The quantity k(s) = |ei s |, 
appearing in Eq. (2), is the curvature (or the inverse radius of curvature) at 
a point s, measuring the deviation of the curve from linearity. With these 
definitions for the unit tangent and unit normal vectors one can complete 
an ordered triad of orthonormal vectors at any point on the curve with the 
introduction of the third unit binormal vector e 3 = e x x e 2 . The derivatives 
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of the unit normal and binormal can be written as 

C2s = -«ei + re 3 (3) 

e 3s = -re 2 . (4) 

Eq. (3) introduces torsion, r = e^e^ x e lss /|e ls | 2 , while Eq. (4) follows 
immediately on account of orthonormality. Just as curvature measures 
deviation of a curve from linearity, torsion measures deviation from 
planarity. A planar curve has zero torsion. 

Eqs. (2-4), constitute the so called Frenet-Serret equations (FSE) 
describing the basic geometry of a curve in 3-D [Stoker, 1969; Struik, 1961]. 
k and r are intrinsic quantities describing the curve. Though the three unit 
vectors can be canonically defined at any point, ambiguity arises at points 
where k vanishes. Thus for a straight line, or along segments of a curve 
that are straight, one does not have a unique choice of e 2 and e 3 . We shall 
assume hereonwards that k does not vanish and is positive. Given the two 
intrinsic parameters k and r, at least in principle, a corresponding curve 
can be constructed, while the actual construction can be much involved. 
This reconstruction problem remains a challenge, except for some simple 
curves, and has deeper significance in soliton theory [Sym, 1985]. Besides 
the abstract interest in the subject, the problem is also of interest in the 
image processing community [Sapiro, 2001]. 

In the same spirit one can also consider describing the time evolution of the 
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curve. First we note that the FSE Eqs. (2-4) can be written compactly as 

e is = x e; ; % = 1,2,3. (5) 

where the Darboux vector Q = Ke^ + rei, or simply (r, 0, k) in the e« basis. 
Let the dynamics of the curve R(s) be given by 

R t = aei + (3e 2 + 7e 3 , (6) 

where now the subscript t stands for time derivative and ct,/3 and 7 are 
arbitrary functions of s and t. If one assumes the curve to be 
non-stretching, then s, and hence |ej|, are independent of time. Though 
this appears to be a serious restriction at first, a variety of systems still can 
be modeled by such curves, while for most other systems this turns out to 
be a good approximation, and definitely a good starting point to familiarize 
with the dynamics of curves. 

The dynamics of the triad can be described by defining a new vector 
uj = (001,002,^3) such that, similar to Eq. (5), 

e it = nxei; i = 1,2,3. (7) 

Note that Eq. (7) can be identified with the equations of a rotating rigid 
body, whose dynamics is described by the three principal axes. The vector 
00 then plays the role of its angular velocity. Consequently, a moving curve 
can be identified with a rotating rigid body moving along the curve 
[Lakshmanan, 1979]. 
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The compatibility condition (ej) st = (ei) ts ,i = 1,2,3, then yields the 

— * 

following equation satisfied by and u. 

(l t - uj s = uj x (l. (8) 

— * 

These are three equations for the five quantities fl and Q. Written 
explicitly, Eq. (8) reads 

K t - u 3s = -tu 2 ; (9) 

T t ~ UJ ls = KUJ 2 , (10) 

u 2s = kujx - tuj 3 (11) 

Any nontrivial choice for u^s in terms of k and r then leads, through Eq. 
(8), to nonlinear time evolution equations for k and r - the intrinsic 
variables representing a curve. Notice that nonlinearity is inherent in the 
system due to the vector product on the right hand side in Eq. (8). The 
compatibility conditions on R, (R st = R is , H sst = Rj ss , etc.,) lead to 

— * 

equations relating a, (3 and 7 with Vt and uj, which however we shall not be 
making much use of here. R(s, t) can also be seen as a surface in 3-D 
generated by a curve. Evidently this surface is special, being generated by a 
curve and one that is non-stretching. 

We list below some important examples of ubiquitous integrable nonlinear 
evolution equations admitting soliton solutions [Ablowitz and Clarkson, 
[1992]; Lakshmanan and Rajasekar, 2003] that can be identified with 
moving curves, for suitable choices of ul as functions of k and r 
[Lakshmanan and Tamizhmani, 1981; Lamb, 1977]. 
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(i) Nonlinear Schrddinger equation (NLS) [Hasimoto, 1972; Lakshmanan, 
1977; Lakshmanan and Tamizhmani, 1981] : For the choice 

ui = r ; u 2 = -k s ; u 3 = -kt (12) 

K 

Eqs. (9) lead to the NLS equation 

| r-s 

iijj t + ^ ss + ttH 2 ^ = ; iJ = nexpi rds'. (13) 

Z J — oo 

(ii) Modified Korteweg-de Vries equation (MKdV) [Nakayama et al., 1992; 
Lakshmanan and Tamizhmani, 1981]: The modified Korteweg-de Vries 

u t + 6u 2 u s + u sss = ; (14) 

is obtained for the choice 

2 2 
Ui = — ; U 2 = ~K S T ; U 3 = -K ss - y+TK, (15) 

with k = 2u and r = constant. 

(iii) sine- Gordon equation [Lamb, 1977; Lakshmanan and Tamizhmani, 
1981]: The sine-Gordon equation 

9 st = sm9; (16) 

is obtained with the choice 

oui — cosm ; u>2 = — sin-u ; u 3 = (17) 

T T 

with k = u s and r = constant. 
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Besides these, many other evolution equations [Lakshmanan et al., 1998] 
and certain integrable hierarchies have been shown to be related to curve 
dynamics [Nakayama et al., 1992; Goldstein and Petrich, 1991]. As we 
noted above, the nonlinearity in the equations is due to the vector product 
in Eq. (8). They continue to be nonlinear even on relaxing the condition of 
local length preservation. In considering applications to certain physical 
systems this becomes inevitable. 



3 Intrinsic Global Quantities Associated 
with Curves 



In this section we introduce a few global quantities, associated with the 
geometry of the curve, which one often encounters and which have 
important physical relevance in applications. The first of these is the 
integral of curvature, the total curvature 



The geometrical meaning of this integral is more apparent from expressing 
the integrand in terms of the tangent vector ei, namely \dei\ = k ds. Since 
ei is a unit vector in 3-D, it can be seen as living on a unit sphere S 2 . For a 
given curve R(s) then, the tangent vector traces a curve, the tangent 
indicatrix, on the surface of this unit sphere [Stoker, 1969]. The integrand 
in Eq. (18) is indeed the length of this curve. For a closed curve in 3-D we 




(18) 
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have a closed indicatrix on the unit sphere (see Fig. 2), and 



N = 



k(s) ds. 



(19) 



However, it should be remembered that the converse is not always true, i.e., 
a closed tangent indicatrix does not always correspond to a closed curve in 
3-D. Through a simple argument it can be shown that, for a closed curve in 
3-D, the total curvature N > 2tc. This is the Frenchel's theorem, which is 
due to the fact that for a closed curve in 3-D, the corresponding curve on 
the unit sphere cannot entirely lie in one hemisphere [Kamien, 2002]. Thus 
the total curvature (or the length of the closed circle traced out on the unit 
sphere) has to be greater that the circumference of the great circle which is 
indeed 2n. 

The second integral we shall consider, and one that is frequently 
encountered, is the 'energy' integral 



We call this the energy integral for this is the first relevant term for any 
typical physical system modeled by a curve. This is understandable, since 
any change in shape of the curve requires energy, whose first nontrivial 
term has to be quadratic. The simplest example is an elastic band, for 
which this is just the potential due to the Hooke's law. Further, if one 
imagines s to represent time and e x as the position vector of a particle, 
then the integrand in Eq. (20) is just its kinetic energy (apart from a 
factor) and E is the action. Since e x moves on a unit sphere, this is just the 




(20) 



11 



action integral of a particle constrained to move on the surface of a sphere. 
E turns out to be the Hamiltonian for the classical continuous Heisenberg 
model for spin systems in 1-D, of which more will be said in section 5. 
It may be naturally tempting at this stage to see the geometrical relevance 
of the other integral, namely that of torsion r 



M 



[ L t ds. (21) 

J 



Analogous to Eq. (18), this can be seen as the length of the binormal 
indicatrix. However, M has a better geometrical interpretation which shows 
up in several practical situations. The answer can be found in Eq. (5) and 

— * 

noting that r is the component of Q along the tangent direction. Thus 

torsion is the local phase acquired by the normal and the binormal as we 

move along the curve. Eq. (21) gives the total phase acquired by the 

(ei — e 2 ) plane as we move through the total length of the curve, compared 

to the same plane in a non-rotating basis. Further, parameterizing the 

tangent vector in terms of spherical polar coordinates, (i.e., 

ei = (sin# cos0, sin^ sin 0, cos#)), we find for a closed tangent indicatrix 

M — jr ds = j> cos 9 d(f). (22) 

Using Stoke 's theorem, the integral above can be rewritten as 

M = <j) cos 9 d(f) = 2n — sin 9d9d(p. (23) 

The last integral in the above equation is just the area enclosed by the 
tangent indicatrix on the unit sphere. Thus, except for an additive 2ti, M 
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gives the area enclosed by a closed tangent indicatrix. We shall encounter 
the same in section 6 when we discuss applications to geometric phases in 
classical optics. In the next few sections we discuss some physical systems 
that can be modeled using space curves, where significant results have been 
reported. 

4 Vortex Filaments 

Dynamics of a fluid is described by its velocity field V. Vorticity is defined 
as the curl of velocity W(= V x V). At low velocities vortices are generally 
observed when a fluid is subject to rotation, the most common example 
being the kitchen sink vortex, while at high velocities, rotation is inherent 
in the fluid. They are known to occur in various forms such as lines, rings, 
etc., [Batchelor, 1967]. In the following we shall confine ourselves only to 
line vortices. Vortices in ordinary fluids have fascinated thinkers over 
several centuries dating back to the masters in Rome. Some early detailed 
observations on fluid vortices can be found in Learnado da Vinci's writings. 
Besides their common occurrence, the geometrical features associated are 
shared by several other systems depicting vortex behavior. These include 
large scale systems such as twisters, magnetic flux tubes in solar corona 
[Spruit, 1981], systems at microscopic scale such as Bose-Einstein 
condensates [Bret in et al., 2003], superfluid helium [Schwarz, 1998] and 
magnetic media. Besides the intrinsic interest in them understanding their 



13 



dynamics has been realized to be the key to fluid turbulence [Frisch, 1996]. 
As is well known, motion of incompressible fluids is described by the Euler 
equations [Saffman, 1992]. For such a fluid V.V = 0, then it follows from 
the Helmholtz theorem that the velocity field at a point X in the fluid is 
derived entirely from its curl - the vorticity W, 



For filament vortices, the vorticity W is along the tangent direction ei, 
W = Cei. We recall two of Helmholtz laws of vortex motion [Saffman, 
1992]: i) Vortex lines move with the fluid, i.e., if the line vortex is described 
by the position vector R, then V(R) = Rt and ii) the strength of the 
vortex tube remains constant, i.e., C is a constant. To a first 
approximation, known as the localized induction approximation, wherein 
one assumes a) thinness of the core radius compared to the mean radius of 
curvature and b) that the vortices are non-stretching [Saffman, 1992], the 
dynamics of these vortices are governed, with a rescaled time, by the 
localized induction equation 



Here, s is the distance along the vortex line. The equation describes the 
motion of the vortex induced by its own geometry. If R is a straight line, 
then evidently the dynamics is trivial. 

Mapping the vortex line on a space curve discussed in section 2, and 




(24) 



R^ — R s x R ; 



'SS > 



Rs-Rs — 1- 



(25) 
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interpreting R s as the unit tangent vector ei, Eq. (25) can be identified as 

e lt = ei x e lss ; ei.ei = 1. (26) 

Now making use of FSE Eq. (5), one can check that Eq. (26) can be 
rewritten as 

e lt = K s e 3 - nre 2 (27) 
One can complete the trihedral evolution as 

e 2 < = -K a ei + (— - r 2 )e 3 , (28) 

K 

e 3 t = -(— - r 2 )e 2 - nre 2 . (29) 

K 

As discussed in section 2, the compatibility conditions e st = e ts leads to the 
evolution equation for the curve that can be written in terms of the 
intrinsic quantities k and r, [Lakshmanan et al., 1976]. 

K t = -(kt) s - K S T, 

r t = [(k ss /k) - r 2 ] s + KK S . 

Further, it was shown by Hasimoto [Hasimoto, 1972] that these equations 
can be combined to a single complex equation - the nonlinear Schrodinger 
equation (NLS) 

^ + ^,.+ 2^ = 0, (30) 

if one defines ip = /texpf/f^ r ds']. The NLS equation being an integrable 
nonlinear evolution equation [Zakharov and Shabat, 1971] (ie., possessing 
an infinite number of independent integrals of motion [Lakshmanan and 
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Rajasekar, 2003; Ablowitz and Clarkson, [1992]]) with soliton solutions, 
similar solutions were conjectured in line vortices also. These were 
subsequently observed experimentally [Hopfinger and Brownand, 1982]. 
The integrable nature of Eq. (25) also guarantees stability of these solitonic 
excitations. 

The form of the vortex filament corresponding to a one-soliton solution of 
the NLS is given in Fig. 3. Explicitly, the expression for the position vector 
R is given by 

R(s, t) = js — yutanh(z/£), — /isech(z/£) cost], (31) 

— yUsech(z/£) sin^j, 

where £ = s — 2\t, r\ = As + {y 2 — \ 2 )t and ji = 2z//(z/ 2 + A 2 ), involving two 
constant parameters v and A. The condition that the dynamics be length 
preserving is a valid approximation for small times. Prominent examples 
where this condition has to be necessarily relaxed include superfluids and 
certain weakly nonlinear oscillations such as in chemical media [Schwarz, 
1998; Cross and Hohenberg, 1993]. Chemical oscillations are governed by 
the complex Ginzburg-Landau equation, 

A t = A(l + ib)VA- (l + ic)\A\ 2 A (32) 

wherein numerical investigations have shown instabilities arising due to 
stretching and bending of the vortex lines [Aranson and Bishop, 1993]. 
Another prominent example where dynamics of vortex filaments play an 
important role is in scroll waves encountered in heart muscles. Scroll waves 
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known to form in the ventricle of the heart in response to certain 
stimulations (Fig. 4). Its geometry can be described by a space curve, the 
filament of the scroll, about which the scroll rotates, with a frequency 
larger than the rate of heart beat. Instabilities in the scroll wave can result 
in a break up into smaller independent scrolls leading to ventricular 
fibrillation which often proves fatal [Winfree, 1987]. 



5 Spin Systems 

Magnetism is a consequence of the quantum nature of the spin angular 
momentum associated with electrons. Specifically, the fundamental 
Hamiltonian responsible for spin-spin interaction that leads to spontaneous 
magnetization in magnetic materials, which is the Heisenberg interaction, is 
a consequence of Pauli's exclusion principle [Kittel, I960]. If 
Sj(= (Si, Sf, Sf)) were to denote the spin vector at a site denoted by '%' , 
then the isotropic Heisenberg Hamiltonian is given by 

H = -JJ2 Si.S,-; |S,| 2 = 1, (33) 

<i,j> 

supplemented by the usual spin commutation relations [Sf, Sj] = 5ij6 abc Sf, 
where Sij and e abc are the usual Kronecker delta and Levi-Civita tensors, 
respectively. Here, the bracket in the sum denotes nearest neighbor. Due to 
the tendency of any system to attain a minimum energy configuration, from 
Eq. (33) we note that the adjacent spins tend to align parallel when J > 0, 
leading to ferromagnetism and anti-parallel for J < 0, giving rise to 
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anti-ferromagnetism. Particular cases of the system lead to other 
interesting models such as the XF-model, n-Potts model and the Ising 
model [Chaikin and Lubensky, 1970]. 

Simple as it may appear, the system, Eq. (33), has posed itself as one of the 
challenging problems in physics. Though further complications to Eq. (33) 
may arise in the form of anisotropy, external magnetic field, dissipation, 
etc., the real difficulty turns out to be due to the dimension and structure 
of the lattice, besides the nonlinearity in their dynamic equations. 
From the Hamiltonian, Eq. (33), using the commutation relations, one 
obtains the dynamical equations for the spins as 

S? t = J^e abc St S], (34) 

3 

where the sum runs over all nearest neighbors of the site i. Though the 
system is discrete and quantum in nature, for low energies it proves 
sufficient to consider the classical continuum version of the model. In this 
case the dynamics of the spin field is governed by a special case of the 
Landau-Lifshitz equation (LLE) given, after appropriate scaling, by 
[Landau and Lifshitz, 1935; Lakshmanan et al., 1976] 

S, = SxV 2 S; |S| 2 = 1. (35) 

For quasi-one dimensional magnets, the evolution equation can be written 

as 

S, = SxS ra ; |S| 2 = 1. (36) 
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Associating now a space curve with the continuum spin chain and 
identifying the spin vector S(x, t) with the unit tangent vector ei and the 
spatial coordinate x with the arc length of the curve s, the evolution 
equation becomes 



Now making use of the earlier identification in section 4 (see Eqs. (27 - 
30)), the one dimensional continuum LLE is seen to be equivalent to the 
integrable NLS, which in turn implies the integrability of LLE in 1-D 
[Lakshmanan, 1977]. 

A spin one-soliton solution is shown in Fig. 5, for which the components are 



where 77, yu. and £ are as defined in section 4 under Eq. (31). 
Besides 1-D ferromagnets, curve dynamics has also been used in studying 
low lying excitations in 1-D anti-ferromagnets, where it is related to a 
special case of static 2-D ferromagnets [Balakrishnan et al., 1993]. Soliton 
solutions in 1-D spin chains are well known. However, the glossary of results 
available in the literature for spin systems in higher spatial dimensions are 
hardly substantial, though several particular cases have been effectively 
dealt with [Kosevich et al., 1990; Bar'yakhtar et al., 1994; Lakshmanan 
et al., 1998]. Notable among these, in the present context, is that of 



e u = ei x e lss ; 




(37) 




(38) 



/zsech(z/£)(z/tanh(z/£) sin r\ — A cos 77), 
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vortex/anti- vortex rings in ferromagnets [Cooper, 1999], and particular 
radially and axisymmetric solutions for the Heisenberg Hamiltonian 
systems [Lakshmanan et al., 1979; Porsezian and Lakshmanan, 1991]. 

6 Geometric Phase in Classical Optics 

The concept of geometric phase in quantum mechanics has gained much 
significance following the work of Berry [Berry, 1984; Shapere and Wilczek, 
1989]. When a quantum system is subjected to a time dependent 
Hamiltonian, the state of the system acquires an extra phase, in addition to 
the dynamic phase, depending on the geometry of the time dependent 
parameter in the Hamiltonian. The simplest example is a charged particle 
in the presence of a time dependent magnetic field vector. Here, the state 
of the particle acquires an extra phase depending on the geometry of the 
magnetic field. Though the concept has achieved much of its significance in 
quantum theory, there are classical examples that do exhibit similar 
features. 

Let us consider the propagation of polarized light along a thin twisted 
optical fiber (see Fig. 6). Electromagnetic waves carry inherently three 
vectors that describe them, namely the electric field E, the magnetic field 
B and the propagation vector k. In free space the three vectors are 
orthogonal, but acquire a component along the other vectors while traveling 
in a dipolar medium. Nevertheless, one can still consider them to be 
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orthogonal to a good approximation, as we shall do here. If one considers 
the fiber to be a curve described by the vector R(s), the direction of the 
propagation vector k is simply the tangent vector e x . The E, B vectors 
acquire a phase depending on the geometry of the fiber, which has been 
experimentally observed [Haldane, 1986; Berry, 1987; Shapere and Wilczek, 
1989; Kugler and Shtrikman, 1988; Chiao and Wu, 1986]. 
Unlike the normal and binormal vectors of section 2 which were canonically 
determined purely from the geometry of the curve itself, the E, B fields 
present a frame that is physically observable. Evidently, they can still be 
expressed in the Frenet-Serret basis formed by e^. An expression for the 
geometric phase acquired can be obtained from the wave equation. 
If R(s) were to describe the center of the fiber, then any point on the fiber 
can be uniquely expressed as 



where r and 9 are polar coordinates described on the circular cross-section 
at s. This holds good for fibers where the radius of curvature at any point 
is much larger than the radius of the cross-section. Defining 
i 2 = cos#e 2 + sin#e 3 and i 3 = ei x i 2 , we find 



It is convenient to transform to orthogonal coordinates (s, r, 0), where 
(f) = 9 + J s t ds' . In terms of these coordinates the gradient operator is 



X = R(s) + r cos 9e 2 + r sin 9e 3 , 



(39) 



dX. = (1 — kt cos 9)dsei + dri 2 + r(d9 + rrfs)i 3 . 



(40) 
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given by 



V = 



1 — KT cos 6 ds 



ei d 



. d '13 d 

dr r dd> 



(41) 



In these transformed coordinates the propagation of the wave, described by 
the wave equation for a field ip, V 2 ip = —Eip, is just that along a straight 
fiber to second order in k [Kugler and Shtrikman, 1988]. Thus the only 
difference in propagation, between a straight fiber and a curved one, is the 
phase factor — 6 = f s r ds'. This is indeed the phase acquired by the 
polarization vector traveling along a curved fiber. 

Though the discussion here was confined to plane electromagnetic waves, 
the result can be generalized to transverse acoustic modes along a curved 
path also, as the conclusion merely follows from the wave equation. Further 
similar geometric phases have also been identified in a dynamic context, 
particularly in continuous ferromagnetic and anti-ferromagnetic spin chains 
[Balakrishnan et al., 1990; Balakrishnan et al., 1993]. 



The last few decades have seen a variety of mathematical methods 
developed and applied to systems that are of interest specifically in life 
sciences. Several filamental structures most common in nature, namely 
biological macromolecules such as proteins, bacterial fibers and DNA, 
tendrils in creeping plants and polymers exhibit complex physical 
phenomena due to the high order of nonlinearity in their governing 



7 Soft Condensed Matter 
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equations [Goriely et al., 2001]. The elementary physical model for most of 
these systems is a thin elastic rod acted upon by forces and moments. 
Besides these, as in the previous sections, we shall also assume a) the 
cross-section to be thin compared to the mean radius of curvature, and b) 
the rod is inextensible and unshearable. Let us identify a local basis along 
the axis of the filament (ei,i 2 ,i 3 ), as defined in the previous section. Recall 
that the last two vectors differ from the normal and binormal by an angle 9. 
The description of the rod dynamics is simplified by averaging out the 
equations for the force F and moment M over the cross-section. After 
appropriate scaling, the model is described by the Kirchoff equations 
[Antman, 1995] 

F ss = ( ei )«, (42) 
M s + e 1 xF = i 1 x {i^a + i 2 x (i 2 ) tt , (43) 

M = a/ciei + « 2 i 2 + bK 3 i 3 . (44) 

Here b = I3/I2 is the ratio of the principal moments of inertia of the 
cross-section and a = J/(2J 2 (1 + a)), where J is a factor depending on the 
particular shape of the cross section and a is the Poisson ratio. For a rod of 
uniform circular cross section, I 2 = I3 = J/2. Eq. (44) defines the moment 
M in terms of the components of the twist vector 
(ki, k 2 , k 3 ) = (r + d9/ds, /tsin^, k cos $) . The model has been used in 
analyzing instabilities arising due to change in curvature, leading to 
spontaneous changes in its conformation. An elastic coil when stretched 
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and untwisted completely and slowly released can go into various possible 
conformations, induced by variation in curvature, such as looping [Goriely 
and Tabor, 1998a], orientation reversal, buckling, etc. Typical examples 
include telephone coils and tendrils in creeping plants [Goriely and Tabor, 
1998b]. 

As we mentioned in section 3, the expression for E, Eq. (20), describes the 
free energy of a flexible string. A more realistic model involves, to 
quadratic order, contributions due to twist and possible external fields, such 
as viscosity: 

E = J (A(k - k ) 2 + Ct 2 ) ds- J A ds. (45) 

Here A and C represent coefficients of bend and twist respectively. k 
represents any possible ground state configuration other than a straight 
rod. This is understandable in case of most polymers or proteins and DNA, 
which all have a nontrivial natural configuration. The second integral 
imposes the constraint of local length preservation. The dynamic equations 
for the elastica are obtained by equating the force — 5E/5K to the viscous 
drag due to the surrounding fluid, proportional to R t . Note that the second 
derivative in time of R does not appear in the force equation, since the first 
derivative is the more dominant term near equilibrium. The model has 
been used in analyzing dynamics of polymers in fluids, leading to geometry 
induced phenomenon such as untwisting and folding [Goldstein and Langer, 
1995; Goldstein et al, 1998]. 
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8 Extension to Surfaces 



The success of the models based on curve evolution prompts us to consider 
systems that can be modeled on surfaces and hyper-surfaces in higher 
dimension. Indeed there has been considerable progress both in the variety 
of physical systems covered in these models and on developing the required 
mathematical methods to deal with these systems [Bishop et al., 1984; 
Pelce, 1986]. In this section we discuss how the ideas we developed for 
curves can be extended to surfaces in 3-D. Let the surface be described by a 
position vector P(x, y) [Struik, 1961; Stoker, 1969]. Just as in the case of 
curves, this allows one to define three vectors (see Fig. 7), namely, 
u = dP/dx, v = dP/dy and a third vector orthogonal to these two, 
N = uxv/|uxv|. Note that unlike the case of curves, the three vectors 
are not orthonormal, in general. While N is orthogonal to both u and v by 
construction, the latter two are neither orthogonal to each other nor 
normalized. For the same reason the normal vector N becomes ambiguous 
when u and v are parallel. We shall consider only smooth surfaces, wherein 
these ambiguities do not occur. The intrinsic variables describing the 
surface are obtained from the first and second fundamental forms, 



respectively. The variation of the three vectors with respect to the two 
parameters is given by the Gauss- Weingarten equations, which again 



{dPf = Edx 2 + 2Fdxdy + Gdy 2 , 



(46) 



-dP.dN = Ldx 2 + IMdxdy + Ndy 2 , 



(47) 
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involve the coefficients of the second fundamental form and derivatives of 
the first fundamental form. The coefficients are not completely independent 
and are related by the Gauss-Codazzi-Mainardi conditions P xxy = P xy x and 
P = P 

1 yyx 1 yxy 

The role of intrinsic curvature for such a surface is played by the Gaussian 
curvature, defined as 

T ^ LN-M 2 
= -EG^- (48) 

For a flat surface, evidently dN vanishes and the Gaussian curvature K 

trivially vanishes. Solitons in surface theory as specific surfaces were known 

to geometers centuries before the formal theory of soliton developed. 

Indeed if one considers a surface of constant negative curvature (say 

K — — 1) -a pseudo spherical surface, the angle 9 between the two vectors u 

and v satisfies the sine-Gordon equation 

9 xy = sin 6. (49) 

The sine-Gordon equation is an integrable model and its role in several 
physical systems was recognized only in the middle of the twentieth century 
[Chaikin and Lubensky, 1970; Jackson, 1990; Lakshmanan and Rajasekar, 
2003]. 

The difference between a smooth surface and one generated by a moving 
curve must be noted. While a smooth surface can be imagined as one 
generated by a moving curve, the converse is not always true. However the 
surface generated by a moving non-stretching curve has the advantage that 
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a triad of unit vectors can be defined at every point, which makes it easier 
to deal with. Besides, the underlying curve being described by two intrinsic 
variables k and r, the surface can be immediately described by a complex 
function such as ip in Eq. (30). 

Consider hence a surface generated by sliding a curve of the type discussed 
in section 2 [Lakshmanan et al., 1998]. Such a surface is essentially 
described by a position vector with two parameters, of the type described 
in the previous paragraph. However this surface is special, as at every point 
on the surface we can define a triad of unit orthonormal vectors 
ei(x, y), i — 1, 2, 3. If one considers the time evolution of this surface, these 
vectors are then described by three sets of equations 

(ei)a. = (I x a ; (ei) y = f x a ; (e;) t = u x a. (50) 

— * — * 

Here T(= (71, 72, 73)) plays a role similar to Q of Eq. (5) and uo of Eq. (7). 
Consequently, along with Eq. (8), we have two more compatibility 
conditions 

Q y - f t = f x u ; r x — Q y = QxF. (51) 

Written in terms of the intrinsic variables k and r, these evidently lead to 
coupled nonlinear partial differential equations in two spatial dimensions 
and one time dimension (i.e., (2 + 1) dimensions). Several (2+l)-D 
generalizations of the Heisenberg equation (33) have been mapped on such 
moving surfaces and further they have been related to integrable (2+l)-D 
nonlinear evolution equations by the above procedure. Some examples are 

27 



as follows [Myrzakulov et al., 1998]. 

(i) Myrzakulov I equation [Myrzakulov, 1987]: 
This equation reads 

S t = (S x S y + uS) x , (52) 
u x = — S.S X x S y . (53) 

Identifying again the spin vector S(x,y,z) with the unit tangent vector 
ei(x,y,z) and following a procedure similar to section 5, we can make the 
following identifications in Eqs. (50). 

72 = -, (54) 

K 

Ul = 1m. - rd^Ty ; uj 2 = -K y ; uj 3 = -nd~ l T y (55) 

K 

where d~ x stands for integration with respect to x. Analogous to the NLS 
equation (30) for the 1-D Landau-Lifshitz equation (36), one can construct 
a complex NLS type equation - the Zakharov equation for a complex 
function q = (/t/2)9 a : 1 (— r) 

iq t (x,y,t) = q xy + Vq ; V x = 2\q\ 2 y . (56) 

(ii) Ishimori II Equation[lshimon, 1984]: 
The equation has the form 

Sf — S X (S XX Syy) + U X Sy + UyS X , (57) 

U XX + Uyy = 2S.S;r X Sy. (58) 
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The procedure adapted for (i) in this case leads to 

^xx ~T ^yy TUJ3 ^2x /rn\ 

= — ; Wl = ^— (59) 

^2 = u x ^ 2 ~ k x + 7 3?/ + 7i7 2 (60) 

^3 = -KT + KU y + M x 7 3 - 7 2j/ + 7273- (61) 

An analogous complex equation can be constructed if one defines p = a/b, 
where 

a = ^ 2 + 7 2 2 + 73 2 - 2*7 2 ] 1/2 (62) 
b = ^ (r - £ + ^-^-^ N (63) 

x v 2 + 72 2 + 7 | - 2/t 72 y v ; 

p satisfies the Davey-Stewartson equation II, 

ipt + Pxx ~ Pyy - 2p<fr = 0, (64) 
<Pxx + (Pyy + (\p\ 2 )xx ~ i\V?)yy = 0. (65) 

Similar connection can be identified between Ishimori-I equation and the 
Davey-Stewartson equation, which are known to possess localized solutions 
called dromions. 

Though examples presented here pertained to generalized versions of spin 
systems, other physical scenarios, such as hydrodynamics, 
magnetohydrodynamics and dynamics of shell membranes, have also been 
effectively investigated using geometry of moving surfaces [Rogers and 
Schiff , 2003] . Other alternative procedures have also been suggested to 
relate nonlinear integrable evolution equations to moving surfaces [Sym, 
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1985; Nakayama and Wadati, 1993; Konopelchenko, 1996]. Notable among 
them is the one introduced by Konopelchenko [Konopelchenko, 1996]. Here, 
one considers surfaces induced in three dimensional space via solutions of 
two dimensional linear problems, while their integrable dynamics is 
obtained by the (2+1) dimensional nonlinear evolution equations obtained 
as compatibility conditions of the linear problems. Nevertheless, the subject 
of relating surfaces to higher dimensional nonlinear partial differential 
equations, and the connection to the larger problem of integrability, 
remains a open and challenging task. 

9 Conclusion 

We have attempted in this tutorial review to illustrate the subject of the 
nonlinear dynamics underlying the curve motion and discussed how this 
arises as a natural model in studying various physically interesting systems. 
That the evolution equations in many cases closely relate to integrable 
nonlinear evolution equations with soliton solutions adds to the interest in 
the subject. This has been illustrated by considering several systems 
wherein curve geometry plays a crucial role. While we have been largely 
concerned with the simpler case of non-stretching curves, the condition can 
be relaxed and generalized to include a wider class of systems often 
encountered. Such studies can also relate the dynamics to interesting 
spatio-temporal patterns. The extension of curve dynamics to surfaces as a 



30 



method to study higher dimensional nonlinear partial differential equations, 
and their integrability, has been discussed as a challenging future direction 
of the ideas well realized in the case of curves. 
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Figure 1: A curve in 3-D represented by the position vector R and 
parametrized by the arc length s. An ordered triad of unit vectors ej,i = 
1, 2, 3 can be canonically defined at any point on the curve. 

Figure 2: A closed curve on a unit sphere traced out by the tangent vector 
corresponding to a closed curve R(s) in 3-D. 

Figure 3: A one-soliton filament vortex in a fluid, corresponding to the one 
soliton solution ip = 2z/sech(z/(s — 2Xt)) expi(As + {y 2 — \ 2 )t). 

Figure 4: A scroll wave. The scroll can be described as a spiral wave around 
a line filament 

Figure 5: A propagating one soliton solution in a one dimensional Heisenberg 
ferromagnetic spin chain. The spins rotate about the z axis. The soliton 
profile describes the spin component in the x — y plane (1 — S 3 ). 

Figure 6: Propagation of polarized light along a curved optical fiber. The 
E, B fields acquire a geometric phase depending on the geometry of the fiber, 
given by j r ds. 

Figure 7: A surface is described by the position vector P(x,y). A pair of 
tangent vectors u, v and a unit normal N can be defined at every point. 



38 



Figure 1 



Figure 2 




Figure 5 



Figure 6 



Figure 




